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R. W. Helmes, 1 M. Sindel, 1 L. Borda, 1 ' 2 and J. von Delft 1 

1 Physics Department, Arnold Sommerfeld Center for Theoretical Physics, 
and Center for NanoScience, Ludwig-Maximilians-Universitdt Miinchen, 80333 Miinchen, Germany 
2 Research Group "Theory of Condensed Matter" of the Hungarian Academy of 
Sciences and Theoretical Physics Department, TU Budapest, H-1521, Hungary 
(Dated: February 12, 2005) 

Recent experiments measuring the emission of exciton recombination in a self-organized single 
quantum dot (QD) have revealed that novel effects occur when the wetting layer surrounding the 
QD becomes filled with electrons, because the resulting Fermi sea can hybridize with the local 
electron levels on the dot. Motivated by these experiments, we study an extended Anderson model, 
which describes a local conduction band level coupled to a Fermi sea, but also includes a local 
valence band level. We are interested, in particular, on how many-body correlations resulting from 
the presence of the Fermi sea affect the absorption and emission spectra. Using Wilson's numerical 
renormalization group method, we calculate the zero-temperature absorption (emission) spectrum 
of a QD which starts from (ends up in) a strongly correlated Kondo ground state. We predict two 
features: Firstly, we find that the spectrum shows a power law divergence close to the threshold, 
with an exponent that can be understood by analogy to the well-known X-ray edge absorption 
problem. Secondly, the threshold energy ojq - below which no photon is absorbed (above which no 
photon is emitted) - shows a marked, monotonic shift as a function of the exciton binding energy 



PACS numbers: 73.21. La, 78.55. Cr, 78.67.Hc 

I. INTRODUCTION 

Recent optical experiments^ 2 using self-assembled 
InAs quantum dots (QDs), embedded in GaAs, showed 
that it is feasible to measure the absorption and emis- 
sion spectrum of a single QD. In absorption spectrum 
measurements photons are absorbed inside the QD by 
electron-hole pair (exciton) excitation. In emission spec- 
trum measurements, on the other hand, an exciton 
created by laser excitation recombines inside the QD, 
whereby a photon is emitted which is measured. 

Due to spatial confinement, the QD possesses a charg- 
ing energy and a discrete energy level structure, which 
can be rigidly shifted with respect to the Fermi energy 
Ep by varying an external gate voltage V s . Therefore 
V g allows for an experimental control of the number of 
electrons in the QD, which in turn determines the en- 
ergy of the absorbed and emitted photons. Indeed, the 
optical data reveal a distinct V^-dependence and justify 
the assumption of a discrete energy level structure of the 
QDi. 

In the experimental set-up, depicted in Fig.fi] the InAs 
QDs are surrounded by an InAs mono-layer, called 'wet- 
ting layer' (WL), like islands in an ocean. Above a certain 
value of V g , the conduction band of delocalized states of 
this WL begins to be filled, forming a two-dimensional 
Fermi sea of delocalized electrons, i. e. a two-dimensional 
electron gas (2DEG). The 2DEG hybridizes with local- 
ized states of the QD, leading to anomalous emission 
spectra which could not be explained by only considering 
the discrete level structure of the QD 1 . 

Motivated by these experiments, we investigate here 
the optical properties of a QD coupled to a Fermi sea, at 
temperatures sufficiently small that Kondo correlations 



can occur (T=0). The Kondo effect in a QD has already 
been detected in transport experiments^, where it leads 
to an enhanced linear conductance. So far the Kondo 
effect in QDs has been studied almost exclusively in rela- 
tion to transport properties. The experiments of Refs. 1,2 
open the exciting possibility to study the Kondo effect in 
optical experiments. 

In optics, the effect of Kondo correlations on QDs has 
to the best of our knowledge been discussed theoretically 
only with respect to non-linear and shake-up processes in 
a QD 5 £. The influence of disorder in heavy-fermion sys- 
tems on the x-ray-photoemission has been studied e. g. 
m Rcf. 7 In this paper we investigate the absorption and 
emission spectra of a QD. We are especially interested 
in optical transitions (examples are shown in Fig. 0] be- 
low) for which the QD starts in or ends up in a strongly 
correlated Kondo ground state, and will investigate how 
the Kondo correlations affect the observed line shapes. 
In Ref£ the emission spectrum in the Kondo regime has 
already been studied, however with methods which only 
produce qualitative results. 

The paper is organized as follows: In Section^ we ex- 
tend the standard Anderson models by including a local 
valence band level (henceforth called v-level) containing 
the holes. In contrast to Refsi^ we consider only one 
local conduction band level (henceforth called c-level), 
to simplify the calculations. In Section II I II we explain 
how Wilson's numerical renormalization group (NRG) 
methodic can be adapted to calculate the absorption and 
emission spectrum of the QD. In Section llVI we present 
the results of our calculations and predict two rather dra- 
matic new features. Firstly, the absorption and emission 
spectra show a tremendous increase in peak height as the 
exciton binding energy t/ cxc is increased. In fact, the ab- 
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back contact 

FIG. 1: Right hand side: experimental setup used in Refs4*- 
(Picture: courtesy of the group of K. Karrai.) [bottom to top: 
GaAs substrate (2000 nm), highly doped GaAs back contact 
(20 nm, zero-point of x-axis), GaAs tunnel barrier (20 nm), 
InAs mono- layer, forming the wetting layer, together with the 
QDs, GaAs layer (30 nm), AlAs/GaAs tunnel barrier (~ 100 
nm), GaAs (4nm), NiCr top gate]. The gate voltage V s , ap- 
plied between the back contact and the top gate, drives no 
current through the system, since the contacts are separated 
by two tunnel barriers. Left hand side: position-dependence 
(in z-direction) and energy-dependence of the lower conduc- 
tion band edge of the layered structure for two different gate 
voltages Vg and V g b ; the different band gaps of each mate- 
rial result in jumps at the corresponding interfaces. Since 
the InAs band gap is smaller than that of GaAs, there is a 
dip in the band gap at the GaAs-InAs interface, resulting in 
QDs with localized conduction and valence band states. The 
number of localized electrons trapped in the QD can be con- 
trolled by Vg, which shifts the energy levels with respect to 
the Fermi energy Ep (set by the back contact). Inset: holes 
can be trapped as well due to the bump of the upper band 
edge of the valence band at the position of the QDs. Ir- 
radiation by laser light excites electron-hole pairs (excitons) 
inside the GaAs layer, which migrate and become trapped in 
the InAs QDs. Finally they recombine by emitting photons, 
whose emission spectrum is detected. 



sorption spectrum shows a power law divergence at the 
threshold energy ujq, in close analogy to the well-known 
X-ray edge absorption problemii. Exploiting analogies 
to the latter, we propose and numerically verify an ana- 
lytical expression for the exponent that governs this di- 
vergence, in terms of the absorption (emission)-induced 
change in the average occupation of the c-level. Secondly, 
the threshold energy below which no photon is absorbed 
or above which no photon is emitted, respectively, say uio, 
shows a marked, monotonic shift as a function of C/ oxc ; we 
give a qualitative explanation of this behaviour by con- 
sidering the interplay of various relevant energy scales. 
Conclusions are given in Section El 
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FIG. 2: Model of a semiconductor QD, consisting of one c- 
level and one v- level, with energies e c and e v , respectively. 
The Coulomb repulsion of two electrons in the c-level has the 
strength U c ■ The coupling between the c-level and the 2DEG 
is parametrized by the tunnelling matrix element V. Crucial 
for the model is the Coulomb attraction between holes in the 
v-level and electrons in the c-level, which has a strength f7 e xc- 
The excitation of electrons from the v-level to the c-level (by 
photon absorption) and the relaxation of electrons from the 
c-level to the v-level (by photon emission) is considered as a 
perturbation of strength 7. 



II. MODEL 

The experimental setup used in Refsi^, which inspired 
our analysis, is depicted in Fig. ^ (see Fig. caption for de- 
tails). To model this system, we consider an Anderson- 
like model 9 for a QD, with localized conduction and va- 
lence band levels, coupled to a band of delocalized con- 
duction electrons stemming from the WL. Our model is 
similar in spirit, if not in detail, to that proposed in 
RefsAiS. It consists of six terms, illustrated in Fig. [21 

H = Ho + H pcl -t, (1) 

where 

Hq = H c + H v + WtJcxc + Hwh + Hc-WL- (2) 

We consider one c-level with energy e c and one v-level 
with energy e v , originating from the conduction or va- 
lence band of the InAs QD, respectively. Note that e v 
is smaller than e c by the order of the band gap; since 
this difference is at least two orders of magnitude larger 
than all other relevant energy scales, its precise value is 
not important, except for setting the overall scale for the 
threshold for absorption or emission processes. 

Since one c-level is sufficient to produce the effects of 
present interest, we will, in contrast to the experimen- 
tal situation realized in Refsi^, disregard further local 
levels to simplify the calculations (the experimental sit- 
uation realized in Refsi^ will be considered in a future 
publication) . 
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The c-level and the v-level are described by H c and 
7i v , respectively, 

cr 

H v = y^eynyo- + U v (l - rtytX 1 - n v ;), (3) 



where 



= rt, 



and 



Here the Fermi 



operators cj. and vj. create a spin-cr electron in the c- 
level or in the v-level, respectively. The parameters U c 
and U v are Coulomb repulsion energies which have to be 
paid if the c-level is occupied by two electrons or if the 
v-level is empty, respectively. Since states with two holes 
are very highly excited states independent of the value of 
U v (due to the band gap), the actual value of U v has no 
influence on the results. The term 



U, 



exc' L C(7 



(l - n va >) 



(4) 



accounts for the exciton binding energy: the Coulomb 
attraction between each electron in the c-level and each 
hole in the v-level lowers the energy of the system by 

U exc • 

The 2DEG formed in the WL is described by 



^WL — E £kl k Jk<r, 



(5) 



where the Fermi operator it creates a delocalized spin-a 
electron with wave vector k. The hybridization between 
the c-level and the 2DEG is described by 



-WL 



k.a 



+ cth 



(6) 



where the tunneling matrix clement V is assumed to be 
real and energy- independent. The hybridization between 
the c-level and the 2DEG is henceforth parametrized by 
r = TrppV 2 , where pf is the density of states (DOS) 
of the 2DEG at Ep; we assume a flat and normal- 
ized DOS with bandwidth D. Since in the considered 
experiments^, the mass of the (heavy) holes is signif- 
icantly larger than the mass of the electrons, we neglect 
the hybridization between the v-level and the 2DEG. 
The last part of the Hamiltonian, 



H 



pert 



E 

k.a 



j k a k e- iuJkt clv a 



(7) 



describing the excitation (first term) and the annihila- 
tion (second term) of excitons in the QD by photon ab- 
sorption or photon emission, respectively, is considered 
as a perturbation of the system. Here a k iflV) destroys 
(creates) a photon of the laser field with wave vector k, 
where the laser has the frequency Lo k — c\k\. The cou- 

pling is given by 7& = e (Jwjk/ZeoV) 2 e k • D, with the 



elementary charge e, the dielectric constant eo, the quan- 
tization volume V, the orientation of the laser field e k , 
and the dipole moment of the QD transition D. Here we 
assume 7 to be independent of k, since we consider the 
laser to be approximately monochromatic. Note that for 
simplicity, we neglect in the present study terms of the 

form J2ka "fk (l\ a v !j + vXho) j describing photon-induced 
transitions between the 2DEG and the v-level. Such tran- 
sitions will lead to Fano-type effects, which we choose not 
to consider here, but will be the subject of future work. 
Treating 7i pcr t perturbatively is valid as long as the op- 
tical line width |7| 2 A where A is the density of states of 
the photon field, is small compared to the Kondo tem- 
perature Tk (defined below), the smallest energy scale in 
our studies: \^\ 2 A < T K . 

In the following considerations, the quantized nature 
of the photon field will not play any role in our consider- 
ations; to calculate emission and absorption line shapes, 
all that we shall be concerned with are the matrix el- 
ements of the operators cj.iv + v' a c a . For simplicity of 
notation, we shall therefore henceforth write the pertur- 
bation term simply as 



Wpcrt =1^2 ( C ^ W,T + V t C °-) 



(8) 



For the scenario of a local spinfull level coupled to a 
Fermi sea, the Kondo effect occurs if the temperature 
T < Tk and the average occupancy of the local level is 
roughly one, i. e. in our case (n c ) = J2 a (n ca ) ~ 1, known 
as the 'local moment regime' (LMR). Here Tk is given 
by 



T K 



(C/ C r/2) 1/2 e ^+u c )/2TU % 



(9) 



see Refii£. If T < Tk, Tk is the only relevant energy scale 
in the problem. The Kondo effect introduces a quasi- 
particle peak, the Kondo resonance, at the Fermi energy 
Ep in the local density of states (LDOS) A c (ui), 



A c (uj) 



E 

{.cr 



(f|4|G) 6(w-(Ei-E & )) 



(Ef Eq)) 



(10) 



see Fig. 13 Here |G) and |f) are eigenstates of Hq with 
energy Eq and Ei, respectively, where |G) is the ground 
state. The LDOS A c (uj) of the c-level is a well-known 
function, which was calculated with the NRG, e. g., by 
Costi et al^, and has been studied frequently since. 

In transport experiments at T < Tk, the Kondo effect 
causes the 'zero bias anomaly', an enhanced conductance 
due to the quasi particle peak at Ep ■ Here we will inves- 
tigate how the Kondo effect affects the absorption and 
emission spectrum^. 

Fig. Ola) and Fig. 0^b) show examples of absorption 
and emission processes to be studied in this paper. For 
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FIG. 3: The normalized— local density of states A c (lo) of the 
c-level in the Kondo regime, with e c = —U c /2. The Kondo 
effect results in a resonance at the Fermi energy E-p. There 
are side peaks of the singly (doubly) occupied local level at 
to = =f[/ c /2 of a level width 2r. 



both examples the QD is tuned such that the c-level is 
initially singly occupied, (n G ) = 1, i. e. in the LMR and 
therefore gives rise to a strongly correlated Kondo state 
for T <T K . 

In the absorption process, Pig. IHa), a photon excites 
an electron from the v-level into the c-level. Due to the 
exciton binding energy, the c-level is 'pulled down' by 
the value of U exc . Thus the occupation of the c-level in 
the final state can have any value between one and two, 
depending on the value of U oxc relative to the charging 
energy U c of the c-level. If the final occupation is not in 
the LMR, the Kondo-state is lost. 

In the emission process, Fig. 0Jb), an electron from 
the c-level recombines with a hole in the v-level, thereby 
emitting a photon. In contrast to the absorption process, 
here the occupation of the c-level decreases since the ex- 
citon binding energy is lost in the final state. Again the 
Kondo state is lost if the final occupation is not in the 
LMR. 



III. 



METHOD 



The absorption and emission spectra can be calculated 
using Fermi's Golden Rule for the transition rate out of 
an initial state |i), which is proportional to 



(E t - Ei)) , (11) 



where |i) and the possible final states |f) are eigenstates 
of Ho, cf. Eq. CJ, with energy E\ and Ef, respectively. 

No analytical method is known to calculate both the 
eigenenergies of Ho and all matrix elements (i|7i pert |f} 
exactly. Here we calculate them with Wilson's NRG 
methodic, a numerically essentially exact method 17 . 




photon 




FIG. 4: QD (cf. Fig. |2J tuned such that the c-level is ini- 
tially (left hand side) singly occupied, (a) Photon absorption 
process, inducing a transition from a state with no hole in the 
v-level and a singly occupied c-level (Kondo state) to a state 
with a v-level hole and a doubly occupied c-level (non-Kondo 
state). As indicated, the occupation of the c-level in the final 
state is determined by the value of e c — U cyLC + U c relative to 
E-p. (b) Photon emission process, inducing a transition be- 
tween a state with a v-level hole and a singly occupied c-level 
(Kondo state) to a state without v-level hole and an empty 
c-level. 



A. Block structure of Hamiltonian 

Since Ho commutes with h vai the number of holes in 
the v-level is conserved. Thus it is convenient to write 
the unperturbed Hamiltonian Ho in the basis | x I')c+wl® 
|4 r ) v , where |\E , ) c +wl denotes a product state of the c- 
level and the 2DEG in the WL, and \^) v denotes a state 
of the v-level. In this particular basis the unperturbed 
Hamiltonian Ho reads 



|0)v 
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H v -\ 



















V o 











Ho 



where the Hamiltonians 

H v o = Hc-wl + Hwl + H c 

7ivT = Hc-wl + Hwl + H c 
H v ^i = Hc-wl + Hwl + H c 



(12) 



(J 
(J 

2e v (13) 



act only on states | v I')c+wl- Since we have not included 
a magnetic field in our model, H v -\ = H v [. 
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Absorption (A) and emission (E) processes (see Fig.0} 
involve transitions between different blocks of Eq. (|12|) : 

A : |G) = |G) c+ wl ® |TI)v |f> = |f)c+WL 8) |a) v 

E : |g> = |g)c+WL ® |CT} V -» |f) = |f)c+WL ® |TI)v, 

(14) 

where |G) is the ground state of Ho, |G) c+ wl the corre- 
sponding ground state of 7i V T| an d |g) c +WL is the ground 
state of H va , with a =|, j. For absorption, which is gov- 
erned by cj.iv, |f) is a state of the block 7i vcr . For emis- 
sion, which is governed by v\c a , |f) is a state of the block 

To calculate the absorption spectrum, cf. Fig. Ufa), 
we insert |G) for |i) in Eq. (JTTJ. Then ac(w) gives the 
probability per unit time for the transition from |G) to 
any final state |f) of Tt VIJ [containing one hole], equiva- 
lent to the probability per unit time that a photon with 
frequency ui is absorbed, which is the desired absorption 
spectrum qig(w), divided by \j\ 2 . The actual value of 
7 is not important, since it does not affect the shape of 
the absorption function, but only its height. The same 
argument applies to the emission spectrum, cf. Fig.^b). 
Here, one needs to insert |g) for |i) in Eq. (JTTJ. 

To employ the NRG to calculate a; (a;) via Eq. (JTTJ, 
one has to overcome a technical problem. The NRG is 
a numerical iterative procedure, where the energy spec- 
trum is truncated in each iteration [besides the first few 
iterations]. In standard NRG implementations, transi- 
tions from or to highly excited states can only be calcu- 
lated qualitatively rather than quantitatively. In our case 
we need to compute transitions to or from states of the 
blocks H v „, see Eq. ifHEf) . which are highly excited since 
they are separated by the order of the band gap from 
states of H v ti, see Section |H] We solve this problem by 
keeping the same number of states for the blocks Ti. VC r 
and 7i v |x in each NRG iteration, which is in principle 
the same as running two NRG iterations for both blocks 
at the same time. This approach is similar to the one 
used by Costi et al^£, who studied a problem analogous 
to ours. 



B. Limiting case of vanishing exction binding 
energy (£/ cxc = 0) 

To check the accuracy of the modified NRG method, 
we begin by considering the limiting case of vanishing 
exciton binding energy, U exc — 0. We will show that for 
this particular case the absorption and emission spectra 
are related to the local spectral function. 

For U cxc = the v-level is decoupled from the c-level 
and the 2DEG, see Eq. 10}. When decomposing the states 
in the same way as above, I*} = |$) c+ wl® \^)v, the total 
energy can be written as a sum, E = £ c+ wl + -Ev Thus, 
using Eqs. and i|14|) . the absorption and emission 



spectrum can be written as 

a G (uj) = 27r^ | c+WL (f|4|G) c+W L| 2 • 

f ,<T 

x 8 (u> - (.E f)C+WL - Bg.c+wl) -Awj , 
a g (w) = 2tt^ | c+WL (f|c (T |g) c+WL | 2 • 

x (E { , c +wL - Sg,c+wl) +Au?j . 

(15) 

Here Aw = E{ :V — Eq, v = — e v represents a constant shift. 

To compare the LDOS with the absorption and emis- 
sion spectrum, we divide it as A c (lo) = A£(u>) + A~(lj), 
with 

At(u>) = |c+WL(f|4|G)c+WL| 2 

x s(uj - (£f, c+W L - £gx+wl)) for io > 0, 

ATM = ^ |c+WL(f|c a -|G) c+W L| 2 

x s(uj + (£f, c +wL - -Eg.c+Wl)) for lo <0. 

(16) 

Since the operator c\ does not change the state of the VB, 
the sum in Eq. (|16|) runs only over states |f) of WvtJ,- 

To compare Eqs. I|15[) with Eqs. I|ltj|) . note that for 
U exc = the blocks of the Hamiltonian l|12|) are degener- 
ate (aside from a constant shift), TL va = T^vU: an d thus 
G) c+ wl = |g)c+WL- Therefore 

A: o-gM = 2ttA+(cj - Aw), 

E: a s (uj) =2-kA-(-u- Aw). (17) 

Thus, for C/ cxc = 0, we can calculate the absorption and 
emission spectra in two different ways: firstly, with the 
modified NRG procedure and secondly, via Eq. (|T7Jl with 
A c (lo) obtained from the NRG as well. We find an ex- 
cellent agreement between both approaches, which serves 
as a consistency check that the modified NRG works as 
intended. 



IV. RESULTS 

In Section Till Bl we showed that for U cxc = the ab- 
sorption or emission spectra are related to the LDOS. 
Starting from this well-understood limiting case, let us 
now study how the absorption and emission spectra be- 
have upon increasing U cxc . We use the modified NRG 
procedure, described in Section ITTT1 to calculate the ab- 
sorption and emission spectra Q!i(w) from Eq. lllj l. The 
results are shown in Fig. |SJ We see two striking be- 
haviors: Firstly, there is a tremendous increase in peak 
height for both the absorption and emission spectra. In 
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FIG. 5: Absorption and emission spectra for different values of t/ e xc (e c /U c = — 0.5,U V /U C = 1.0 and T/U c = 0.15). An increase 
in t/exc results in an increase in the slope of the divergence at the threshold-energy loo and in a monotonic shift of u>q of the 
absorption and emission spectra. The right panel shows the divergence of the emission spectrum at the threshold energy u>o 
normalized to T K [T K /U C = 0.020, extracted from Eq. ©]. 



fact, we find that the spectra diverge at the threshold en- 
ergy loq , the energy below which no photon is absorbed or 
emitted, respectively, in close analogy to the well-known 
X-ray edge absorption problem. Secondly, the threshold 
energy u>o shows a marked, monotonic shift as a function 
of the exciton binding energy U exc . 



A. Exponent of the power- law divergence 

Let us first study the divergence of the spectral peak at 
threshold. For any U exc 7^ 0, we find a power-law diver- 
gence for both the absorption and the emission spectra—, 
for energies u> near the threshold energy loq: 
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a(co) ~ , u> -> uj . (18) 

\UJ - UJqJ 

Examples of this behavior are shown in Fig. H3 where 
the absorption spectrum is plotted for several different 
values of U cxc on a double logarithmic plot leading to 
nice straight lines for energies (lj — luq) < Tk, i. e. 
in the regime where Kondo correlations can build up. 
The slope of such a line yields the exponent (3. Re- 
markably, we find that the exponent so determined de- 
pends only on the change An in the local occupation, 
An = ± ({n c )f — (n c )i) ('+' for absorption, '— ' for emis- 
sion) , to be called 'screening charge', and obeys the fol- 
lowing relatioiiSS (who's origin will be discussed below): 



FIG. 6: Asymptotic behavior (lj — > cjq) of the shifted absorp- 
tion spectra normalized to Tk [Tk/U c = 0.10, extracted from 
Eq. ©]. For energies (uj — luq) < Tk, in the regime where 
Kondo correlation build up, we find the power-law behavior 
as predicted in Eq. 1181 . The exponent (3 increases as (7 CX c is 
increased. The lower bound of ckg(w) is set by the number 
of NRG iterations (here: (lo — loo) ~ 10~ 3 Tk). The asymp- 
totic behavior of the emission spectra (not shown here) looks 
identical to that of the absorption spectra. 

^An-M. (19) 

Since for u> — > luq the relevant transitions in the case of 
absorption and emission are |G) — ► |g) and |g) — > |G), re- 
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FIG. 7: The exponent /3 of the power-law divergence ex- 
tracted from the NRG results for different values of e c (sym- 
bols) coincides very well with the formula for /3 given by Eq. 
1191 (solid line), indicating that (5 is fully determined by An. 
Here An has been varied between and ~ 0.8 by varying U cxc 
between 0U C and 1U C in steps of 0.1J7 C - 



spectively, An is the same for both types of transitions, 
implying the same exponent (3 for both absorption and 
emission for a given choice of parameters. In particular, 
for the absorption spectra whose asymptotic behavior is 
shown in Fig. [U we have (for uj — > u> ) (n c )i = (n. c )c 
and (n c )f = (n c ) g where (n c ) g and {n c )c denote the av- 
erage occupation of the states |g) and |G), respectively, 
see Eq. (|14|) . At J7 CX c = 0, we have An = 0, since there 
|G) c +wl = |g)c+WL, see Section As U cxc increase, 
(n c ) g and thus An also increases, since the Coulomb at- 
traction between the hole and the electrons in the c-level 
pulls down the c-level to an effective value e c = e c — U cxc 
[note that |g) is an eigenstate of Tt va , whereas |G) is an 
eigenstate of Ti. v ]i and thus independent of U exc \. 

We have extracted the exponent (3 for several different 
values of e c . For each value of e c we have varied An 
between and ~ 0.8 by varying C/ e xc between and U c . 
The results are shown in Fig. [7| We find a very good 
agreement between the results extracted from the NRG 
and the universal behavior predicted in Eq. H19|) : all data 
points nicely collapse onto the curve predicted by Eq. 

The numerical results presented in Fig. [S] should thus 
be interpreted in the following way: for U exc — we have 
An = and thus (3 = 0, which gives a finite height of the 
absorption and the emission spectrum at the threshold 
[in fact the height is 2ir times the height of the corre- 
sponding LDOS, see Eq. I)17p]. As soon as we choose 
values of {7 CX c > 0, we find (3 > 0, leading to an infi- 
nite height of the absorption and the emission spectral 
peaks. Of course, the infinite peak height is not resolved 
by our numerical data, for which a(uj — luq) is always fi- 
nite. However, with increasing U cxc the exponent (3 also 



increases, resulting in a steeper slope of the peak at the 
threshold, which leads to a higher peak in the numerical 
results. 

An explanation for the universal behavior given by Eq. 
(|19|l can be given by studying the analogy between the 
physics presented in this paper and the well-known X- 
ray edge absorption problem. A result analogous to Eq. 
(|18|l was found by Schotte and SchotteSI, where the ab- 
sorption spectrum was studied for the X-ray edge prob- 
lem. [In RefiSi all results are presented for the absorp- 
tion spectrum. However, by rewriting Eq. (7) of Rcf. 21 
for emission, their results can be applied to the emission 
spectrum as well. Keeping that in mind, we will focus 
only on the absorption spectrum in the following, but 
the argumentation can easily be applied to the emission 
spectrum as well.] In Ref— , Fermi-liquid arguments re- 
lating phase shifts and local screening charges are used 
to derive an expression for the exponent (3, namely 



(3=1-^1 



(20) 



where N a is the 'effective' number of spin-cr electrons [not 
necessarily an integer] which flow away from the local 
level in the absorption process. Eq. (|20|l is known as 
"Hopfield's Rule of Thumb"—. We can use this result 
to analyze our absorption spectra, too, since the system 
behaves like a Fermi liquid for T = 0. Thus arguments 
based on the relation between phase shifts and screening 
charges do apply. In experiments, we expect to find the 
behavior J5DJ| for A\j\ 2 ,T < to - lu < T K , where A\j\ 2 
is the optical line width. 

To see that Eqs. ijUlj) and i(2Tj|) are equivalent, we will 
now analyze the absorption process |G) — > |g) [relevant 
absorption process at threshold] and count the charges 
N a . It is helpful to consider an example for the process 
G) — > |g), shown in Fig. [51 where the initial state |G) is 
the strongly correlated Kondo ground state with a singly 
occupied c-level. The state |G) is a coherent superpo- 
sition of states with different occupation of the c-level, 
where the contribution of the state with empty c-level is 
small but finite [depicted in Fig. IHtb)] . If the operator 
c\v a [the part of 7i pcr t corresponding to absorption] is 
applied to |G), this contribution results in a state with 
one hole, a singly occupied c-level and one extra spin-cr 
electron in the conduction band, illustrated in blue in Fig. 
|S{a). The latter subsequently flows away from the QD, 
making a contribution of +1 to N a . This contribution 
to the final state c\v a \G) also is a part of the state |g), 
which likewise has contributions from states with empty, 
singly and doubly occupied c-level [Fig.|S{c)]. The weight 
of the contribution with singly occupied c-level to |g) de- 
pends on U cxc : the Coulomb attraction of the hole in 
|g) pulls down the c-level to the effective value e c result- 
ing in an increase of the average occupation (n c ) g of the 
c-level by An compared to (n c )c see above. As J7 CX c is 
increased, the charge An [which screens the Coulomb po- 
tential of the hole] increases and thus the relative weight 
of the contribution to |g) with doubly occupied c-level 
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small, but finite weight relative weight depends on U 



FIG. 8: (a) Illustration of an example for the absorption process |G) — * |g), the relevant absorption process for energies uj 
close to the threshold wo for the case e c = —[4/2, for which the average occupations of the v- and c-levels are (ti v }g = 2 and 
{n c )G = 1. Figs, (b) and (c) schematically depict the initial and final states |G) and |g), respectively. Both |G) and |g) are 
coherent superpositions of very many different components, whose c-level can be either empty, singly or doubly occupied; for 
|G), the v-level is doubly occupied, and the components with empty and doubly occupied c-level have a very small but finite 
weight. For |g), the v-level is singly occupied, and since the electron-hole attraction lowers the energy of each c-level electron by 
[4xc, the weight of the components with singly or doubly occupied c-level depends on [4 XC . m the absorbtion process depicted 
in (a) and described by the matrix element (g|c^« CT |G), a photon causes the promotion of a spin-er electron from the v-level to 
the c-level (middle panel). This leads to a transition from a |G)-component with empty c-level and an extra electron in the 
conduction band (which subsequently flows away from the QD, making a contribution of +1 to A4), to a state with one v-hole 
and a singly occupied c-level, which in turn is a component of |g) [see (c)]. The Coulomb attraction in the state |g) between 
the v-hole and the c-level electrons pulls down the c-level from e c to an effective value e c = e c — [4xc, resulting in an increase 
of the average c-level occupation (n c ) g by An compared to (n c )a- The screening charge An = (n c ) g — {n c )a flows towards the 
QD, thus making a contribution An/2 to both A4 and A4 (with a = {j, f } for a = {j, J.}). 



increases, too, whereas the relative weight of the state 
with singly occupied c-level decreases. The screening 
charge An flows towards the QD, making equal contri- 
butions -An/2 to both N a and N s . (with a = {|, f} 
for a = {Tjl})- [Another possibility for a transition 
form |G) to |g) [not depicted in Fig. Ufa)] starts from 
a component of |G) with a singly occupied c-level and 
ends up in a contribution of |g) with doubly occupied 
c-level. One obtains the same results for N a and N s if 
one argues that one unit of charge with spin a has to 
leave the doubly occupied c-level and the charge An has 
to flow into the c-level to reach the average occupation 
(n c ) g = (n c )c + An.] Collecting all contributions to N a 
and N s , we find N a = 1 - An/2 and N s = -An/2, 
which, when inserted into (|2U[) . yields Eq. i|19|) . 



A similar argument has been used iniS^, where the 
local spectral function of the Anderson was studied. 



B. Behavior of the threshold energy luq 

Let us now consider the second effect observed in Fig. 
03 the monotonic shift of the threshold energy luq. The 
threshold energy for both absorption and emission is 
given by uo = E g — Eq, where Ho\G) = Eq\G) and 
Ho | g) = -E g |g), as explained in Section ITTT1 The shift in 
wo can be understood by considering a mean-field esti- 
mate of the relevant energies Eq and E g : 

i\ \ 2 

-E G - 2e v + e c (n c ) G + [/ c ^-n c ^) G , 

/l \ 2 

E g ~ e v + e c (n c ) g + U c (^-n c J -U exc (n c ) g . (21) 

Here a correlation energy of the order of 7k has been ne- 
glected. The average occupations (n c )c and (n c ) g can be 
calculated by NRG. Eq. (|21|l allows for a rough estimate 
of the threshold energy ujq: 

UIq = E g — Eq 
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FIG. 9: Behavior of the threshold energy as a function of 
C4xc- Upper panel: average occupation (n c ) g of the c- level 
for the state |g), see Eq. I|14|l . Three distinct regimes can be 
identified: empty orbital (I), LMR (II) and full orbital (III) 
regime, where the c-level is empty, singly or doubly occupied, 
respectively. Since the Coulomb attraction between a hole in 
the v-level and the electrons in the c-level 'pulls down' the 
c-level [e c — e c ~ U axc ], an increase in U BXC bears the same 
effect for (n c ) g as a decrease in e c . Lower panel: threshold 
energy uio versus t/ GX c extracted from the NRG results (solid) 
and obtained from the mean-field estimate (circles), Eq. 122L 
where (n c ) g from the upper panel has been used. 

~ -e v + e c ((n c ) g - {u c )g) 

+-U C {{n c )l - (n c )%) - U cxc (n c ) g . (22) 

The results for ujq shown in Fig. [§| reveal a good agree- 
ment between the threshold energy extracted from the 
absorption and emission spectra calculated with NRG 
(solid line) and the estimation given by Eq. (|22|l . In the 
latter approach (n c )Q and (n c ) g were determined via the 
NRG, see top panel of Fig. [note that (n c )c does not 
depend on U cxc and that (n c )G = ("c)g for U exc =0]. We 
find a linear behavior of ojq as a function of U exc for those 
values of U QXC where (n c ) g stays approximately constant. 
For this purpose three regions of constant occupation can 
be identified, region I {(n c ) g ~ 0), II {(n c ) g ~ 1) and III 
(("c) g ~ 2). As expected by considering the last term 
in Eq. (J22J), we observe the slope of wo(£4xc) to be in 
region I, to be —1 in region II and to be —2 in region III, 
respectively. 

The cross-over regions (dotted lines in Fig. where 
(n c ) g changes between and 1 (I — > II) or between 1 and 
2 (II — > III), on the other hand, show non-trivial behavior 
as a function of uj - In these regions the terms in Eq. © 
compete with each other, which explains the non-linear 
behavior (since in these regions (n c ) g itself is a function 



of U CXC7 too). 

V. CONCLUSIONS 

Motivated by experimental studies of excitons in QDs 
coupled to a wetting layeri^, the aim of this paper was to 
calculate the absorption and emission spectra of a QD in 
the strongly correlated Kondo ground-state. We studied 
an extended Anderson model, including a local valence 
band level and a local conduction band level which is 
coupled to a Fermi-sea (2DEG), see Section [n] For the 
academic limiting case of a vanishing exciton binding en- 
ergy, C/cxc = 0, we could relate the absorption and emis- 
sion spectrum to the well known local density of states of 
the local conduction band, see Section ITTT1 Starting from 
this limiting case, we used the NRG to study the spectra 
for arbitrary values of U cxc . Our main results are summa- 
rized in Fig. [S]which shows two rather dramatic features: 
Firstly, an increase in the slope of the divergence of the 
absorption and emission spectrum as U exc is increased. 
In fact, the spectra show a power-law divergence at the 
threshold energy. Remarkably, the exponent of the diver- 
gence depends only on An, the difference in occupation 
of the local conduction band level between the inital and 
final states for transitions at the threshold. We showed 
that the universal behavior of the exponent can be ex- 
plained by considering the X-ray edge problem, which 
stands in close analogy to the physics presented in this 
paper. Secondly, increasing U cxc produces a marked shift 
of the threshold energy, which can be understood rather 
simply on a mean field level. 

In the present paper we considered, for simplicity, 
a model which contains only a single local conduc- 
tion band level. However, in the present generation of 
experiments^, the wetting layer forms a 2DEG only for 
values of the gate voltage V g for which several local con- 
duction band levels are occupied (not only one, as as- 
sumed in the present paper). Nevertheless, we expect^, 
that future generations of samples could be produced for 
which the assumptions of our model, namely one c-level 
with presence of a 2DEG, are fulfilled. Of course it would 
be very interesting to generalize our considerations to 
more general models, including several local conduction 
band levels. 
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